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I. 


ANALYSIS 


We  consider  the  problem  of  inferring  the  most  probable  location 
of  a target  given  measurements  of  its  photon  emission  collected  from  a 
satellite  equipped  with  a number  of  planar  photon  detectors. 

Specifically,  we  assume  that  the  satellite  carries  nj^  detectors 

whose  respective  orientations  are  given  by  the  angle  of  elevation 

from  the  nadir  and  the  azimuthal  angle  a measured  counter-clockwise 

from  the  satellite's  direction  of  motion.  The  count  rates  from  the 

detectors  are  recorded  at  times  t^,  i = 1,  Z,  . . . , n.^,.  This  provides  a 

total  of  ng  = n^^  . n.^  positive  numbers  designated  by  C^,  i = 1 n^. 

Because  the  C.  include  an  unknown  factor  depending  on  the  strength  of 

the  target  we  replace  then  by  a corresponding  set  of  scaled  quantities 

S.  -having  unit  mean.  Thus, 
i 

S.  = /t(C.  - Cg):  K = ng/Z(C.  - Cg)  (1) 

where  C„  is  the  background  count  rate  (assumed  isotropic). 

The  position  of  the  satellite  is  assumed  to  be  known  at  each  instant 
t..  Thus,  for  a given  target  location  it  is  clearly  possible  to  compute  the 
scaled  quantities  purely  from  the  geometry:  each  will  be  inversely 

proportional  to  the  distance  to  the  target  and  directly  proportional  to  a 
known  function  of  the  angle  subtended  at  the  detector  by  the  target  and 
the  normal  to  the  detector  face.  For  an  ideal  detector  this  function  is 
just  the  cosine,  but  an  empirical  response  function  given  in  tabular  form 
can  equally  well  be  used.  The  actual  location  of  the  target  is  taken  to  be 


the  one  tlial  minimizes  tlic  differ<;nces  lietween  the  computed  and  measured 


I 


1 


I 

J 


S.  . The  computation  takes  the  form  of  a series  of  successive  approxima- 
tions starting  from  an  initial  'guessed'  target  location.  This  initial  location 
is  conveniently  taken  to  be  at  the  point  on  the  surface  of  the  earth  that  is 
viewed  by  the  detector  with  the  highest  count  rate. 

We  describe  the  location  of  the  target  by  giving  its  latitude  X,p.  and 
its  longitude  /i,p.  Then,  in  vector  notation 

S = f(n  (2) 

where  S = js.,  i = 1,  n^j  , F=  |Xp,,  | and  f is  a non-linear  function 

given  by  tht  georoetry.  Equation  (2)  is  approximated  locally  by  the  linear 
equations 


s-f(ro>  ^ c.{FQ)Ar  (3) 

where  is  the  current  approximation  and  AF  the  first  correction,  and 

G(F)  = 9 i/df'  is  a x 2 matrix  of  forward  differences.  For  example 


T2 


(4) 


Equations  (3)  constitute  a set  of  ng  algebraic  equations  in  2 
unknowns  whose  least-squares  solution  forms  the  basis  for  the  iterative 
solution  of  the  non-linear  equations  (2) 
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(5) 


T 

Here  the  subscript  k denotes  the  kth  iteration,  G is  the  matrix  transpose 
of  G and  C ^ is  the  inverse  of  the  covariance  of  the  measurements.  Since 
the  variance  is  equal  to  the  mean  for  a Poisson  process  and  since  the 
count  rates  are  uncorrelated,  C is  a diagonal  matrix  with  on  the  main 

diagonal. 

It  remains  to  determine  the  form  of  in  Eq.  (2)  from  which  the 

elements  of  the  matrix  G can  be  calculated.  Let  r^,  r^,  r^^  denote 
respectively  the  radius  vectors  irom  the  center  of  the  earth  to  the  target, 
the  satellite,  and  the  intersection  of  the  surface  of  the  earth  with  the  line- 
of-sight  from  the  detector.  Then  for  each  estimate  f. 


1.  = K cos 
1 


4 TSD  = K 


(6) 


where  K is  a normalizing  factor  which,  for  any  given  set  of  n^  measure- 
ments, is  determined  as  in  Eq.  (1).  Here  the  cosine  function  is  used,  but, 
as  already  mentioned,  it  can  equally  well  be  replaced  by  an  empirical 
function  in  tabular  form. 

Clearly 

r.j,  = j y^,  z^j  = jcos  cos  cos  X^  sin  Mp,  sinXp  j 
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To  find  Tg,  it  is  convenient  to  express  them  first  in  a primed 

coordinate  system  such  that  r is  along  the  x'  axis  and  the  angular 

o 

velocity  vector  of  the  satellite  along  the  z'  axis.  Then 


r's  = 3«'5.  y'5  z'5!  = (a  hg)  j 1,  0,  0 


= ^'£)'  y'^)'  ^ sinccos  a , 


SD  s in  f s in  O' 


where  a is  the  radius  of  the  earth  and  SD  = jr^  - rgj  is  readily  obtained 
by  squaring  the  vector  relation 


"S  " = "d 


and  solving  the  resulting  quadratic  to  get 


2 2 2 1/2 
SD  = (a  -t-  hg)  cos  € - [a  - (a  (-  hg)  sin  « ] 


Finally,  the  values  in  the  unprimed  system  can  be  obtained  by 


r = T r' 

S S 


'•d  = '^^'d 
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where 


T = 


u cosa -(-  u sina,  v cosa  + v sina,  w cosa+  w sina 
X y ’ X y X y 

- u sinCf  + u cosa,  - V sina  ■(-  v cosa,  -w  sina+  w cosa 
X y X y X y ' 


w 


Here  u,  v,  w are  a triad  of  unit  vectors  given  by  the  satellite  ephemeris 
program  which  specify  the  position  and  velocity  of  the  satellite  in  an 
inertial  coordinate  system,  and  a is  the  corresponding  angle  through 
which  the  earth  has  rotated.  In  particular,  u points  towards  the  satellite, 
V is  along  its  horizontal  component  of  motion  and  w is  perpendicular  to 
its  orbital  plane. 

The  result  given  by  Eq.  (6)  must  be  checked  in- each  case  to  make 
sure  the  target  is  not  over  the  horizon  with  respect  to  the  satellite.  From 
the  geometry  it  is  clear  that  if 

HI  = cos  4 TOS  - sin  (6  - ) 


where  cos  6 = a/ (a  + h.^,)  and  sin  ^ - a/ (a  + hg),  then  HI  < 0 implies 

the  target  is  over  the  horizon.  Similarly  negative  values  for  F in 

Eq.  (6)  indicate  radiation  entering  the  detector  from  the  rear.  In  both  of 

these  cases  the  corresponding  entries  in  the  matrix  are  set  to  zero. 

Here  h_  and  h are  the  heights  of  the  satellite  and  target  above  the 

surface  of  the  earth.  (Unless  input  differently,  h^  is  assumed  to  be 

2 

50  km, above  which  height  there  is  about  1 g/cm  of  air.  ) 
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II.  Description  of  data  input  and  output. 

The  data  are  read  by  subroutines  DATIN,  START  and  EPHM. 
The  latter  subroutines  are  called  from  DATIN  and  read  the  orbital 
element  sets  provided  by  GWC;  DATIN  reads  the  input  cards  that 
specify  the  detector  positions  and  readings.  A complete  data  set 
consists  of  the  following: 


CARD 

COL 

VALUE 

DESCRIPTION 

FORMAT 

1 

1-10 

ERR 

Background  count  rate 

floating  point 

11-20 

EFLAG 

= 0 (or  blank),  means  read 
orbit  element  set 
^ 0,  new  flight  not  specified 

floating  point 

21  -22 

MAX  IT 

number  of  iterations 

integer 

(cards  2 and  3 are  read  when  EFLAG  = 0 and  are  the  2-line 
orbit  element  set  formatted  by  Space  Defense  Center) 


2 

3-7 

ISAT 

satellite  number  integer 

19-20 

lYR 

epoch  year  (last  2 digits)  integer 

21-32 

DAY 

epoch  day  (day  + fraction  floating  point 

of  day  of  year) 

34-43 

XNDOT 

1 st  time  derivative  of  mean  floating  point 
motion  or  Ballistic  coefficient. 

45-52 

XNDDOT 

2nd  time  derivative  of  mean  floating  point 
motion. 

3 

9-16 

XIO 

inclination  (deg) 

18-25 

ASNO 

right  ascension  of  ascending  floating  point 
node  (deg) 
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CARD 

COL 

VALUE 

DESCRIPTION 

FORMAT 

27-33 

EO 

eccentricity 

floating  point 

35-42 

ARGO 

argument  of  perigee  (deg) 

floating  point 

44-5  1 

XMO 

mean  anomaly  (deg) 

floating  point 

53-53 

XNO 

mean  motion  (rev/day) 

floating  point 

4 

1 -5 

KDAY 

6-10 

KMON 

day,  month  and  year  (last 
2 digits)  of  data  acquisition 

integer 

11-15 

KYR 

5* 

1-2 

ND 

number  of  detectors 

integer 

AZ 

EL 

azimuth  and  elevation  of 
detectors,  read  in  pairs 
(deg) 

floating  point 
(8F7.  0) 

6" 

1-2 

NT 

number  of  detector 
readings 

integer 

TIMES 

time(s)  (seconds)  of 
detector  readings 

floating  point 
(8F7.  0) 

7* 

SM 

NT  X ND  detector  readings 
input  so  that  the  first  NT 
readings  are  for  the  first 
detector,  etc. 

floating  point 
(8F10.  0) 

8 

1 

D 

= IHD,  program  calculates 

alphabetic 

initial  estimate  of  target 
location 


= IH  , use  following  values 
as  initial  guess 


T* 
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CARD 

COL 

VALUE 

DESCRIPTION 

FORMAT 

2-11 

TARLOC  (1) 

latitude  (deg)  of  target 

floating  point 

12-21 

TARLOC  (2) 

longitude  (deg)  of  target 

floating  point 

22-31 

TARLOC  (3) 

altitude  (km)  if  different 

floating  point 

from  preset  value 

♦ 

Additional  cards  may  be  necessary  for  these  inputs. 

The  program  ends  when  a blank  card  is  encountered  as  the  first 
card  of  a data  set. 

The  output  of  the  program  at  the  kth  iteration  will  consist  of  the 
current  estimate  of  the  target  location  X._,,  representing  latitude 

f and  longitude  in  degrees  together  with  the  determinant  of  the  covariance 

of  error  matrix 

COV  = >c  [g"^  ( r^)  . . G{ 

^ where  K is  defined  in  Eq,  (1)  and  G,  and  C after  Eq.  (5). 

2 1/2 

Also,  the  residual  R = [H(S^  - f^)  ] is  printed  at  each  iteration. 

Note  that  the  procedure  tries  to  minimize  the  positive  definite  form 

Q(X^  , = E(S.  - f.)^/S. 

giving  the  weighted  sum  of  squares  of  the  residuals.  In  some  cases  the 
I iteration  scheme  may  fail  to  converge  if  the  data  and/or  the  initial  guess 


are  too  poor.  In  such  cases  it  is  helpful  to  have  an  isoresidual  contour  plot 
obtained  from  a tabulation  of  Q ( A U over  a grid  of  values  of  X j* 


Such  a tabulation  is  provided  at  the  end  of  the  printout.  From  it  one  can 
determine  if  any  well  defined  minima  exist,  and  if  so  select  the  most  likely- 
one  as  the  initial  guess. 

III.  Sample  Problem 

Following  are  the  printouts  from  two  sample  runs.  The  first  uses  ideal 
data  to  test  the  program,  and  the  second  uses  simulated  real  data.  Four 
detectors  are  placed  at  azimuth,  elevation  angles  of  (45  , 30  ),  (135  , 30  ), 

(225°,  30°),  (315°,  30°);  the  target  is  placed  at  latitude  30°,  longitude  20°. 

In  the  first  run  count  rates  are  input  corresponding  to  3 satellite  locations 
taken  2 seconds  apart.  The  count  rates  are  proportional  to  the  number 
obtained  by  computing  the  cosine  of  the  angle  between  the  target  and  the 
detector  normal  and  then  dividing  it  by  the  square  of  the  range.  The  back* 
ground  count  rate  is  taken  to  be  zero.  From  the  printout  it  is  clear  that 
the  program  finds  the  target  after  2 iterations. 

In  the  second  run  a set  of  physically  realistic  count  rates  are  selected 

from  Poisson's  distribution  having  the  count  rates  of  the  first  example  as 

means.  This  is  accomplished  approximately  by  replacing  the  count  rate 

c.  bv  c.  t z.yF  where  Z.  is  taken  from  a table  of  normal  deviates 
i ^ i i I I 

with  zero  mean  and  unit  variance.  A background  count  rate  of  10  is  assumed. 
We  see  that  the  program  comes  to  within  1 degree  of  the  target  after  2 


iterations. 


I V! 


APPENDIX 
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IMF  IVAN  A.  OKI  riNCi  I^UK)HA  lOKIKS 


n f I .»b.»ratc>r\  Operations  of  Ihe  AeroHp.ue*  f ..rporalion  is  torjduMtny 
ryj-f  M m«- nt.il  nrul  theoretit  al  j nvr  st  i gat  i ons  nmessary  for  the  evaluation  anrl 
amih'  ation  of  siientitu  ,id\  anc  es  to  new  military  <on<,i  pts  and  systems.  V«*r 
satilitv  anrl  flexibility  base  beer.  d«-veloped  to  a huih  <leiire»-  by  the  latioratorv 
personnel  in  dealing  with  the  many  problems  encountered  -n  the  nation's  rapidK 
de.elopin^’  spai  e and  missih*  systiuns,  Kxpertise  in  the  falest  s^ientifu  devcl 
opments  is  vital  ti>  the  at  c ompl  i shment  of  tasks  relatefl  to  th«*se  problems.  He 
l.ib..rat«>ries  that  <ontribule  to  this  research  are: 

Ae  rophysir  s f .aboratory:  Launch  and  reentry  a<- rodynamic  s,  heat  t rans - 

fer.  reentry  physics,  chemical  kinetics,  structural  mechanics,  fliijht  clyriatnics, 
atmospheric  {Killution,  and  high  • powe  r gas  lasers. 

<!^hemistrv  and  Physics  Laboratory:  Atniosph<*n<  reactions  and  atrnos- 

pl.er.c  oplL  8,  cheroical  reactions  in  polluted  at  n^osphe  re*  s,  chemical  reactions 
of  excited  species  in  rocket  plumes,  chemical  t he  rmodvnann  < s.  plasii.a  and 
laser  - indue  ed  reactions,  laser  chemistry,  propulsion  chemistry,  space  vacuum 
and  radiation  effects  on  materials,  lubrication  and  surface  phenomena,  photo- 
sensitive materials  and  sensors,  high  precisum  lascT  ranging,  and  the  appli- 
cation of  physic  s and  < hetni  st  ry  to  p roblems  c>f  law  enfor<  enu*nt  and  biomedicine , 

kdectronns  K«»searrh  Laboratory:  k'.le<  t romagnetic  theory,  devices,  and 

pr«.pagatic>n  phenomena^  includrng  plasma  eU«c  t romagnetic  s.  quanturri  electronu  s. 
lasers,  and  el<*ct  ro- opti  c s ; communication  sciences,  applied  electronics,  semi 
c «>nduc  ting,  superronduc  tmg.  and  crystal  device  physics,  optical  and  acoustical 
imaging,  atmospheric  pollution:  millimc*t<*r  wave  ancl  far-infrared  technology. 

Materials  Sri«*nc«»s  Laboratory:  Development  of  n»*w  materials,  nictal 

matrix  composites  ancl  new  forms  of  carbon,  tc'st  and  c*.aluatioTi  of  graphite 
and  c e ramu  8 in  reent  ry  , spac  ec  raft  mate  ri  a|  s and  elec  t roni  c component  s i n 
nuclear  weapons  environment,  application  of  fracture  mcmhanii  s to  stress  cor- 
rosion ancl  fatigue*  - indue  ed  fractures  in  structural  r^ietals. 

Spac  c*  Sc  lenc  es  Laboratory:  Atmospheric  and  ionospheric  physics,  radia- 
tion from  the  atmosphere,  dc*nsily  and  composition  of  the*  atmosphere,  .lurorac* 
and  airglow.  magnetospheric  physics,  cosmic  rays,  g«*neration  ancl  propagation 
of  plasma  waves  in  the  magnetosphere;  solar  physics,  stiidic*s  of  solar  magnetic 
fic'Ids.  spate  astroncmiy,  x-rayasf  ronomy;  the  effect.**  of  nuc  lea  r explo.sions, 
magnetic  storms,  and  solar  activity  on  the  earth's  at  nmsphe  re , lonospher-  , and 
magnetosphere,  the  effects  of  optical,  clc*ct  roiTiagnetic , and  particulate  radia- 
tions in  spac  e on  spac  e systems. 

TiiK  akhospac;k  COKPOKATION 
k'.l  Segundo,  ( alifornia 


